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MODULATION OF DYNAMICAL SYSTEMS FOR COMMUNICATION 

PRIORITY INFORMATION 

This application claims priority from provisional application Ser. No. 60/269,052 filed 
5 February 15, 2001. 

BACKGROUND OF THE INVENTION 
Many approaches have been taken to embed continuous-time information bearing 

£3 waveforms onto continuous-time carrier waveforms for the purpose of communication. 

O 

^ 10 Amplitude modulation (AM) and frequency modulation (FM) are common examples of such 
approaches to communication. In both AM and FM, the carrier wave is a sinusoidal wave. 
There is a need in the art for a systematic procedure for constructing continuous-time modulators 
and their corresponding demodulators using carrier waves that are generated by nonlinear 
systems that have periodic, ahnost-periodic, quasi-periodic or chaotic attractors. 
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SUMMARY OF THE INVENTION 

Accordingly, the invention provides a new class of signal modulators and their 
corresponding demodulators. The modulator embeds into a carrier signal an information signal 
by modulating the oscillatory rate of the carrier signal in a manner proportional to the 
20 information signal. The permissible carrier signals are any signals that can be generated by a 
nonlinear dynamical system that has a known exponentially convergent observer. 

According to one aspect of the invention, a method is provided of embedding an 
information signal into a carrier signal that is generated by a non-linear dynamical system by 
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modulating the rate at which the non-linear system evolves in a manner proportional to the 
information signal. The method includes multiplying the information signal by a constant to 
produce a first signal value. The metiiod further includes adding the first signal value and a 
nominal rate of evolution of the dynamical system to generate a second signal value. The 
method further includes providing a feedback path that includes a first and second path, wherein 
the input to the first path is the integration of a multiplication of the second signal value and 
output of the second path, and the input of the second path is output of the first path, such that 
the second patti is a first function that defines the non-linear dynamical system, and providing the 
output of the first path as input to a second function that produces a transmitted signal. 

According to another aspect of the invention, a system for demodulating a transmitted 
modulated carrier signal is provided. The demodulator system includes an observer component 
that receives as input the transmitted signal and a rate estimate and produces an estimate of the 
current state of the transmitter. The demodulator system also includes a rate estimator that 
receives as input the transmitted signal and the estimate of the current state of the transmitter to 
produce the rate estimate. The observer component and rate estimator are interconnected in a 
feedback loop arrangement, and as shown below, this arrangement recovers the information 
signal assuming that the dynamical system possesses a local exponential observability property. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of the modulator system; 

FIG. 2 is a block diagram of the demodulator 20; 

FIG. 3 is a block diagram of the demodulator 30; 

FIG. 4 is a schematic of an exemplary modulator circuit 40; 



FIG. 5 is a schematic of the observer component 50 of the exemplary demodulator circuit 

52; 

FIG. 6 is a schematic of the rate estimator component 60 of the exemplary demodulator 
circuit 52; and 

FIG. 7 is a table illustrating exemplary components of the exemplary modulator and 
exemplary demodulator circuits 40, 50, and 60. 

DETAILED DESCRIPTION OF THE INVENTION 
FIG. 1 is a block diagram of the modulator system 10. An information signal m(t) is 
supplied to a modulating system 10. The modulatmg system 10 consists of a non-linear system 
that has a periodic, quasi-periodic, almost periodic, or chaotic attractor. The information signal, 
m(t), is embedded iti the non-linear system 10 by modulating the rate at which the dynamical 
system evolves. The non-linear dynamically system used to construct the system in 10 is 
denoted using the notation 

x = tyj(x), (1) 
where x is an N-dimensional vector and o)^ is a constant. Modulation is achieved by applying 
the information signal, m(t)^ to the non-linear dynamical system as 

x = {(o,^Pm{t)){{x\ (2) 
where co^ is the nominal rate of evolution and )S is a parameter that characterizes the degree of 
modulation. The effects of modulation on the non-linear system are readily seen by denoting the 
solution to the system without modulation giving in (1) as x^t) ^ Thus, the addition of 



modulation as described by (2) results in the modulated signal that can be expressed in term of 
x„it) as 

x(t) = x,(t + fil[m{r)dt). (3) 

Modulating in this fashion does not modify the attractor of the non-Unear system. It only 
modulates the rate at which the dynamical system evolves along the attractor. 

The signal that is transmitted to the receiver is a scalar function of the state variables of 
the non-linear dynamical system in the transmitter. As shown in FIG. 2, y(t) is denoted as 

yit) = h(t), (4) 
where h(x) is a function that maps the N-dimensional vector, x, to a one-dimensional signal, y. 
The complete modulator system is represented mathematically as 

i = i(0,+fim{t))f(x\ 
y = h(x), 

where x is a N-dimensional signal and j is a one-dunaisional signal, 

FIG. 2 is a block diagram of the demodulator 20. The demodulator 20 must track 
changes in the rate of evolution of a dynamical system. The demodulator 20 consists of two 
fundamental components - an observer component 22 and a rate estimator component 24. The 
observer component 22 reconstructs the state of the demodulator 20, z , given y(t) and the rate 
estimate, th(t) . It is assumed assume that the dynamical system without modulation as given in 
(1) has a known exponentially convergent observer. We define an exponentially convergent 
observer to be a dynamical system 

z = f(z,A(Xo)) (6) 
that has the property that ||z-Xo||<e"'^for some 1>0, where Xq is the solution to the 
unmodulated dynamical system given in equation (1). The rate estimator estimates m(t) given 



the state estimate, z, and the transmitted signal, y. The intercomiection of these components 
comprises the complete demodulator system 20. Given the dynamical system is used in 
modulator system 10, the observer component 24 is given by 

z = (fi>,+)5m)f(z,j;), (6) 
where m is the rate estimate. The rate estimator 24 takes as input the reconstructed state z from 
the observer component 24 and the transmitted signal y(t) and tracks m(t). The low-pass filter 26 
removes any spectral energy known to be absent from the original modulating signal, m(t). The 
low-pass filter 26 is an optional component to the demodulator 20. 

Assuming that the unmodulated dynamical system used in modulator 10 has a known 
exponentially convergent local observer function, the observer component 24 can be modified so 
that it is an exponentially convergent observer of the modulator 10 when m(t) = /«(, , where is 
an unknown constant. If the rate estimator 24 converges to a value of , then the augmented 
observer is assumed to be able to track a time-varying m(t) provided that m(t) varies 
sufficiently slow. 

The design of the rate estimator 24 is based on a technique that is referred to as a 
backwards perturbation expansion. The essential step in this perturbation expansion is to express 
the modulator state, x , as a perturbation expansion about the demodulator state, z , in terms of 
the rate estimate error, e„ =m-mo. By expanding the modulator state x as perturbation about 
the demodulator state z , the resulting e3q)ansion terms depend only on variables local to the 
demodulator 20. The perturbation variables can then be combined in such a way that they force 
the demodulator rate estimate error to zero. If a dynamical system used in a modulator is of the 
form 



(7) 
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and has an exponentially observer of the form 

then a modulator as can be constructed as 

y = h(x). 



(8) 



(9) 



which is similar to the modulator 10 defined in (5) above and a demodulator will be 



constructed as 



y = h{z), 

^=K{y-y)g[^(^H^y 



(10) 
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dz 



which is the same as demodulator 20 as defined in (6) where 0<K<K* for some K*>0, and 
10 r{; •) is a scalar valued function defined by 

/at \ ^ ^ 

(11) 



where ^(0 is scalar function such that sgn(g(a)) = sgn(a) , and 



(12) 



for some y>0. The observer 22 and rate estimator 24 are intercormected in a feedback loop as 



15 shown in FIG. 3. 
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Many systems have been shown to possess an exponentially convergent observer, 
including some chaotic systems. One such system is a chaotic system that is described by the 
Lorenz equations that are defined as 

^i=(o-(^2-^i)) 

= (rXj-x,j^-X2) (13) 

where cr , r, and b are constant parameters. An exponentially convergent observer of the system 
given in (13) when yit) = x^it) , 

z,=cr(z2-Zi) 

z,={ry-yz,-z,) (14) 

FIG. 3 is block diagram of a demodulator 30. Three enhancements have been made to 
the demodulator 30. First, a low-pass filter 36 is added between the rate estimator 32 and 
observer 34 to remove spectral energy in m(t) that is known to be outside the bandwidth of 
m{t) . Second, K is increased beyond the value for which the perturbation expansion analysis 
guarantees stability. Although stabihty is no longer guaranteed, the demodulator 30 may remain 
stable over a range of K>K\ A larger value for K allows mit) to track faster signals. Fmally, 
many of the nonlmearities in ttie demodulator 30 are removed by approximating the rate 
estimator 32 with a linear system. The system that results from the approximation is equivalent 
to a least-squares approach to designing a demodulator. 

To remove the spectral energy in m(t) that is outside the bandwidth of m(t), a filter 36 is 
added to the feedback path. The filtering operation, denoted as <•> , is given by 
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where t/(r,T) is the filter kernel. If the support of this kernel is sufficiently small compared to 



the rate at which m and m vary, tiien 



{^) = -K{ 



— (z)-^, -g -T-i^)-^! [m-m) 
,dz J \dz J I 



S 



(16) 



5 As mentioned above, the gain parameter K had to be smaller than some K* to guarantee 

that the perturbation expansion was bounded. For K<K\ the rate estimate converges 
exponentially and monotonically. Choosing K sUghtly larger than tC affects the demodulator 
30 in two ways. First, since K scales the derivate of m , increasing K also increases the rate 
which m can vary, allowing m to track signals that vary more rapidly. Second, the perturbation 

1 0 term may not remain bounded. Returning to the equation for the perturbation variable. 
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where 



Kz,^x) = 



(dh. . A 



(dh, . , 



(18) 



If K is set to zero in (18), remains bounded. However, the dynamics of demodulator 
15 30 is changed and convergence is no longer guaranteed, demodulator 30 may be stable for a 
range of K>K* . 
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An example, consider the Lorenz based modulation/demodulation system based on the systems 
described in (13) and (14) with K>K* . The term Kr{z,^^) is dropped from (17) and the 
demodulator equation for the Lorenz equation becomes 

h={ci>,+pm){yz2-bz,). (19) 
w = A:(Zi-j)sgn(^,), 
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The modulator equations are 

X2 = (o)^ + fim){rx^ - x^x^ -x^) (20) 

The demodulator 30 has additional nonlinearities added beyond those aheady present in 
the dynamical system. These additional nonlinearities appear in the equation for %^ as given in 

(17). Even when K>K* and Kr(z,%J is removed, a nonUnear equation remains. 

The last term, /3f{z,y), also appears in the observer component 34 of the demodulator 

30. Since f (z, j;) is required by the observer component 34, removing it from the rate estimator 
32 does not reduce the total number of nonlinearities present in the demodulator 30. This term is 

left as it is. The first term, (o)^+j3m) — ^j, is generally nonlinear and does not appear else 

dz 

where in system 30. Approximating this term with a linear, time-invariant system simplifies the 
hardware implementation of the demodulator 30. 
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First, if 0)^ » pm then 

%,^a>,^{z,y)%,-pi{z,y). (20) 
oz 

The varying gain matrix, — {z,y) , is generally nonlinear. However, the differential equation in 

(20) is linear with to and is a time varying linear filter with pf{z,y) as its input Using the 

notation (•) to denote the filtering operation, 

^,^-fi{f(z,y)). (21) 

Replacing this linear time- varying filter with a hne^ time-invariant filter makes the equation for 
consists of only linear components and f (x, j) , the latter of which is already present in the 
demodulator 30. 

The difference between the derivatives of j) and y can be approximated as, 

dh dh 
y-y = i6),+/3m)—(2.)fiz,y)-{a), + ^m-e^))— (z+4e„ +...)f(z+^ie„ +...,y) 

dz dz (22) 

^ej^{z)i{z,y). 

oz 



Filtering y-y with the same filter that appears in (21) results in 

{9-y)^pi^im^,y)y., (23) 

assuming that varies slowly with respect to the time constant of the filter so that can be 
moved outside the filtering operation. Combining (21) and (23), the rate estimator 32 becomes 

^ = Kfi{S^-y)l^{z)i{z,y)) 



«Kp{^{z)fiz,y)^ (m-m). 
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(24) 



From (24), the rate estimator 32 has the form 



(25) 



where a(t) is positive semi-definite function, which suggests m converge to m , 

FIG, 4 is a schematic of a modulator circuit 40. The modulation/demodulation technique 
of the invention is apphcable to any non-hnear system meeting the requirements described in the 
previous section, a particular implementation which uses a set of equations known as the Lorenz 
equations is further described below. The Lorenz equations are a set of non-linear ordinary 
differential equations that have state trajectories that behave chaotically. A signal that is chaotic 
is one that exhibits long-term aperiodic behavior that has a sensitive dependence on initial 
conditions. 

The Lorenz equations discussed above are repeated here for the circuit, hi this circuit 

implementation, however, care must be given to ensure that the signal levels remain within the 

11 2 
operating range of the circuit components. Thus, JSfj =— Xj, ^2 =— and X^=—x^. The 

2 2 9 

rescaled Lorenz equations implemented are 



X,^{o>,^ ^m^X,X,-hX,\ 



(26) 



The implementation of this modulator 40 using multipliers, operational amplifiers and 
capacitors, and resistors is shown in FIG. 4. Solving for the circuit voltage gives 
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It is assumed that the output of the multipliers is the product of its input divided by 10, which is 
typical of multiplier circuits. The nominal rates at which the modulator circuit 40 modulates are 
governed by the capacitors. 

FIG. 5 is a schematic of the observer component 50 of demodulator 52. After solving the 
circuit voltages in the observer portion of the demodulator gives 
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where y = . 

FIG. 6 is a schematic of the rate estimator component 60 of demodulator circuit 52. After 
solving the circuit voltages in the rate estimator 60 portion of the demodulator gives 



m = - 
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where o)^ = — - — and (•) denotes the filtering with a first order low pass filter with a cut-off 

"^309^02 

frequency of o)^^ . The circuit of FIG. 7 utilizes the multiplier that is particul^ to the AD374. 
The AD374 has four inputs. It multipUes the difference between inputs one and two with the 
difference between inputs three and four. The multiplier M3 multiplies y minus the output of 
0P2 with the output of OPl. 
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FIG. 7 is a table illustrating exemplary components of the modulator 40 and demodulator 
circuits 50 and 60. The components values are illustrative of the circuits described in FIGs. 5, 6, 
and 7. Based upon these component values 

=400,000 

y5 = 80,000 (18) 
=80,000. 

However, the modulator and demodulator circuits 40, 50, and 60 may utilize other component 
values. 

The goal of the invention is to describe and analyze a class of signals that can be used in 
the framework of the modulation technique described above, which is referred to as a 
generalized frequency modulation. The invention has various practical application which are 
low power communication system, because the modulation technique is applicable to nonlinear 
dynamical systems, thus the system is not constrained to operate circuit components in their 
linear regime. This potentially reduces the number of circuit components, simplifies the circuit, 
and increases the efficiency. Another is the Spread-Spectrum communication system, which the 
modulation technique of the invention can be applied to chaotic systems, which are naturally 
spread-spectrum signals. Due to sensitive dependence on initial conditions exhibited in chaotic 
systems, chaotic signals are difficult to track without precise knowledge of all of the parameters 
of the chaotic dynamical system, which suggests that chaotic carrier signals may be 
advantageous in the context of private communications. 

Although the invention has been shown and described with respect to several preferred 
embodiments thereof, various changes, omissions and additions to the form md detail thereof, 
may be made therein, without departing from the spirit and scope of the invention. 

What is claimed is: 
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